Stroboscopic observation of quantum many-body dynamics 
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Recent experiments have demonstrated single-site resolved observation of cold atoms in optical 
lattices. Thus, in the future it may be possible to take repeated snapshots of an interacting quan- 
tum many-body system during the course of its evolution. Here we address the impact of the 
resulting quantum (anti-)Zeno physics on the many-body dynamics. We use the time-dependent 
density-matrix renormalization group to obtain the time evolution of the full wave function, which 
is then periodically projected in order to simulate realizations of stroboscopic measurements. For 
the example of a one-dimensional lattice of spinless fermions with nearest-neighbor interactions, we 
find regimes for which many-particle configurations are stabilized or destabilized, depending on the 
interaction strength and the time between observations. 
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Introduction. In the last few years ultracold atoms 
in optical lattices have proven to be a versatile tool for 
studying various quantum many-body phenomena [U |2] . 
Recently, tremendous progress has been achieved by im- 
plementing single-site resolved detection |3l |4] and ad- 
dressing of atoms. Taken to the next, dynamic level, 
one may envision observing the evolution of nonequilib- 
rium quantum many-body states via periodic snapshots 
revealing the position of each single atom. For simpler 
systems, the effect of frequent observations on the decay 
of an unstable state (or on the dynamics of a driven tran- 
sition) has already been discussed and observed, leading 
to the notion of the quantum (anti-)Zeno effect [6-9 . 
Zeno physics has also been seen in cold-atom experi- 
ments with atomic loss channels [10 and was theoreti- 
cally addressed in pJTfM l. Experiments with single-site 
detection, however, would reveal the effect of observa- 
tions on the dynamics of a truly interacting quantum 
many-body system. Here we exploit a numerically effi- 
cient approach to simulating the repeated observation of 
many-particle configurations in interacting lattice mod- 
els. This represents an idealized version of the dynamics 
that may be realized in future experiments. We elabo- 
rate the main features of this "stroboscopic" many-body 
dynamics in the case of a one-dimensional (ID) lattice 
of spin-polarized fermions with nearest-neighbor interac- 
tions. We find a variant of the quantum Zeno effect and 
discuss its tendency to inhibit or accelerate the break- 
up of certain many-particle configurations. In particu- 
lar, the decay rate of these configurations depends in a 
nonmonotonous fashion on the time interval between ob- 
servations. Later, we show that a similar behavior is ex- 
pected for the Fermi-Hubbard and Bose-Hubbard model. 
The discussed features may be seen, e.g., in the expan- 
sion of interacting atomic clouds in a lattice. 

Model. In this paper, we study spin-polarized fermions 



in a ID lattice governed by the Hamiltonian 

n = -jY^ (4q+i + H.c.) + ^ ^ n,n,+i . (1) 

i i 

The first term describes hopping with amplitude J be- 
tween adjacent sites, the second encodes the interaction 
between fermions at neighboring sites, with hi = c-Ci. 
The Hamiltonian displays a dynamical V —V symme- 
try which shows up in expansion experiments [15 . Fol- 
lowing analogous steps as in [15 , we can conclude that 
if both the initial state and the experimentally measured 
quantity O are invariant under both time reversal and tt- 
boost (a translation of all momenta by tt), the observed 
time evolution {0{t)) is identical for repulsive and at- 
tractive interactions of the same strength. The initial 
occupation number states and the n-particle density ob- 
servables in our case fall within the scope of this theo- 
rem. Thus, the only relevant dimensionless parameters 
are \V/J\ and the rescaled time between observations, 
JAt. 

Single particle. We first briefly turn to the single- 
particle case, in which Eq. ([T]) leads to a tight-binding 
band E{k) = —2J cos{k). A particle located initially at a 
single site is in a superposition of all plane wave momenta 
k = — TT . . . TT. After a time t, the probability of detecting 
it at a distance / from the initial site is p(/, t) = J7^^| {2Jt)^ 
where J7 is the Bessel function of the first kind. This is 
shown in Fig. [TJa). The particle moves ballistically, with 
(P) = 2(Jt)^. When the particle is observed repeatedly, 
at intervals At, the ballistic motion turns into diffusion. 
In this case, after m time steps of duration At = t/m, we 
have (P) = 2J^tAt. Thus the motion slows down, and 
in the limit of an infinite observation rate, the particle is 
frozen, which is known as the quantum Zeno effect. 

Simulation. Ideally, each observation is a projective 
measurement in the basis of many-particle configurations 
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Fig. 1. (Color online) Time evolution of the density pro- 
file of fermions expanding in a lattice without observation, 
(a) for a single fermion; (b)-(d) for 13 fermions, initially lo- 
cated at adjacent lattice sites, calculated for several interac- 
tion strengths. 



(occupation number states in real space). To generate 
such measurement outcomes in a numerically efficient 
way, we start by randomly drawing the position of the 
first particle from a distribution given by the one-particle 
density. Afterward, we draw the position of the second 
particle, conditioned on the location of the first one, and 
proceed iteratively (see Supplemental Material [16] for 
details). For the interacting many-body case to be dis- 
cussed now, we use the time-dependent density-matrix 
renormalization group (tDMRG) [ 17H2Q] to calculate the 
joint probabilities and the time evolution between ob- 
servations. We choose a time step of J5t = 0.1 and a 
lattice of typically 115 sites, and we keep up to approx- 
imately 1000 states, at a truncation error of 10~^. For 
the dynamics of noninteracting fermions, an exact for- 
mula can be used (see Supplemental Material [E]). Note 
that tDMRG has been employed recently for dissipative 
dynamics of cold atoms [2T[ |^ . 

Stroboscopic many-body dynamics. We will focus on 
the expansion of an interacting cloud from an initially 
confined state. For 2-species fermions, such an expansion 
has been recently observed in 2D [151 and theoretically 
discussed for ID [23 . We will first briefly address the 
evaporation itself and then discuss qualitatively the re- 
sulting stroboscopic dynamics, with a more refined anal- 
ysis presented further below. Figures [TJb)-(d) show the 
effect of the interaction on the unmeasured time evolu- 
tion of the density profile. For increasing interaction the 
fermions tend to remain localized near their initial posi- 
tions. For interaction strengths \V/J\ > 3 and the times 
shown here, Jt < 16, a more detailed analysis reveals 
that evaporation proceeds via the rare event of a single 
fermion dissociating from the edge of the cloud. The par- 
ticle then moves away ballistically. This evaporation pro- 
cess is hindered by the formation of bound states. This is 
a crucial phenomenon that we will also encounter in the 
context of repeated measurements. For smaller interac- 
tion strengths (|1^/J| < 2), the fermions split gradually 
into a larger and larger number of clusters as time in- 
creases. The parameter regimes in which the model ([T]) 
exhibits diffusive or ballistic transport was addressed us- 
ing tDMRG in Ref. [24^. 
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Fig. 2. (Color online) Specific realizations of the expansion of 
13 fermions with site-resolved detection during the evolution. 
We show the full evolution of density even between observa- 
tions, which collapse the many-body wave function at regular 
time intervals J At [indicated by the dashed lines in panel 
(h)]. Without interaction [panels (a)-(d)], the initial decay 
rate of the configuration increases for larger J At, while for 
large interaction \ V/J\ the decay rate is biggest for finite J At 
[(g), (i)]. For small J At the dynamics becomes independent 
of |y/J|,see (a),(e),(i). 



The effects of stroboscopic observation are shown in 
Fig. [2j for typical realizations of this stochastic pro- 
cess. For noninteracting fermions we find the behavior 
expected from the single-particle case. The spread (and 
thus, the diffusion constant) increases with larger obser- 
vation time intervals J At. For very small J At (strong 
Zeno effect), the motion is diffusive with a small diffusion 
constant that becomes independent of In general, 

it is useful to discuss the initial decay rate of the cluster 
that evaporates via expansion. For the interacting case, 
this decay rate is largest at some finite observation time 
interval J At [Figs. [2|g) and|2|i)], while it is reduced for 
large J At [Figs.[2|h) and|2|l)]. Apparently, at very large 
|y/J|, the initial decay rate may have yet another local 
minimum for intermediate J At, see Fig. [2|j). We con- 
firm this striking nonmonotonous behavior of the initial 
decay rate by simulating 400 realizations for each panel 
shown in Fig. [2] and plot the average number of fermions 
at the central 15 lattice sites as a function of time, in 
Fig. [sj For sufficiently large this number decays 

roughly linearly at a rate that sets the initial decay rate. 
We will see that these features can be mainly attributed 
to a bound state and the two-level dynamics between the 
initial state and the state with a fermion detached from 
the others. 

Doublets and the role of interactions. The effect of 
interactions can be discussed already for the stroboscopic 
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Fig. 3. (Color online) For the expansion process of Fig. 2] the 
average number of fermions remaining at the 15 central lattice 
sites is shown as function of the evolution time Jt for different 
observation time intervals J At. Thin dotted lines correspond 
to expansions without observation, (a) Without interactions, 
the initial decay rate increases monotonically with J At. (b) 
At VjJ — 3.8, the decay rate first increases, then decreases 
with JAt. (c) At VjJ — 9, the decay rate is nonmonotonous 
even for intermediate JAt; compare JAt = 0.3, 0.7, 1.0. Note 
that the lines in (a)-(c) are almost identical for JAt = 0.1. 

dynamics of two fermions. We focus on the decay of a 
doublet, i.e., two fermions sitting at neighboring sites. 

In the quantum Zeno limit, JAt <Cl (or JAt <C |2 J/l/| 
for large see below), only single hopping events occur 
during At. The probability for a fermion hopping left or 
right in this time interval is 2(JAt)^. Thus, the average 
decay time of a doublet is (Jt) = 1/(2 JAt), independent 
of V . This result holds also for clusters of more fermions, 
where the leftmost and rightmost fermion dissolves with 
probability (JAt)^ during At [cf. Figs. [2]; a), [^Je), and 

For larger JAt, the interaction will become impor- 
tant, which gives rise to a bound state (this effect also 
exists for clusters of more particles [25^). Considering 
the basis \l,K) = ^.e:Kjp{iK[j + //2]}c]c]^Jvac) of 
the 2-particle sector with (positive) relative coordinate 

center-of-mass (cm.) coordinate j + //2, and total 
wavenumber K = {ki ^ /c2)mod27r, the action of the 
Hamiltonian ^ is iL\l,K) = \K)^iLK\l)- The first part 
describes a plane wave with wavenumber JT, the second 
the relative motion given by 

Hk\1) = -2Jk [|/ + 1) + (1 - - 1)] + V6i^i\l), (2) 

with JC-dependent hopping amplitude Jk = J cos{K/2). 
A bound state exists \V\ > |2Jx|. It is given by 
li^x) '^'^i{—'2JK/Vy~^\l). We now discuss the de- 
cay of a doublet [see Fig. Qa)] by means of the doublet 
survival probability Pd(^), i-e., the probability of finding 
the doublet intact after time t. Without observations. 



PD{t) = = l,L'\e 



-int\ 



l,L)p, where L and 



L' are cm. coordinates. Thus, in the limit t ^ 00, we 
find Pd(oo) = ^ /q dK\{^ljK\l = 1)1^ Specifically for 
large interaction \V/2J\ > 1, we have 

PD{oo) = l-{2J/Vf + ^{2J/V)\ (3) 

While Pd (00) is determined by the bound state, the evo- 
lution for times Jt < 1 can be approximated by the two- 



doublet Hp^p^-^ 

\/^9"^9^ 9^\/^9^ bound 0.9 \ \ 
state o.8| ^ 1 



0.1 

0, 



V/J=6, JAt=0.6 



0.8 
0.6 
0.4 
0.2 
0.0, 







- - - - 


; V/J 

--..9 




— 6 




••••3.8 




" — 2 






0.5 1.0 1.5 2.0 



Jt 



1 

Pd 
0.6 

0.4 

0.2 

0.0, 



6 7 Jt 





T 




f ^ 


[♦3.8 - 


X 

M « 




♦ 




t ♦ 
114 



0.5 1.0 1.5 2.0 JAt 



Fig. 4. (Color online) (a) Doublet decay level scheme. The 
doublet is separated from the continuum of unbound states by 
an energy gap V. (b) Probability Pd of finding the doublet 
intact after evolution time Jt. Dashed horizontal lines show 
Pd(oo) found in Eq. (|3|. (c) Single trajectory of Pd for a 
time evolution subject to observations, (d) Doublet survival 
probability Pd as function of the observation time interval 
JAt > 0.02 for a fixed total evolution time Jt = 18. Note the 
non-monotonous dependence on JAt for finite interactions. 



level dynamics between |/ = 1) and |/ = 2). This gives 



PD{t) = l 



1 r 

Jo 



dK 



cos2(i^/2) . 



C2 



^m\2iKJt), (4) 



with = [(£^)^ + cos2(Jr/2)]^^^In the strongly inter- 
acting regime we find three regions for the doublet sur- 
vival probability: for times Jt <C ^^=0 probabil- 
ity is independent of the interaction strength, Poit) = 
1 — 2(Jt)^; for times Cx=o ^ < 1 one expects an 
oscillating behavior of Pd {t) given by Eq. Q with a pe- 
riod approximately ^ for |T^/4J| 1; and for Jt ^ 1 
the probability approaches Pd(oo). The full evolution of 
Poit) using exact diagonalization is shown in Fig. pFb). 
Poit) is interaction- independent at times Jt < 
Temporal oscillations in Poit) develop for higher inter- 
action strengths {V/J > 3.5). These oscillations suggest 
that in the presence of stroboscopic observations, illus- 
trated in Fig. [4]^c), the survival probability will depend 
nonmonotonically on the observation time interval. This 
effect is confirmed in Fig.|4]^d). In that figure, the obser- 
vation time interval JAt is varied, while keeping the total 
evolution time constant, Jt = 18 (with a corresponding 
number of observations t/At). The stroboscopic evolu- 
tion is interaction independent for small JAt. For larger 
JAt there is a drastic recovery of Pd in the strong in- 
teracting case, which can show oscillations as a function 
of JAt. This behavior agrees with the one of clusters of 
more fermions, see Figs. [2]and|3j and does not depend in 
detail on the total time Jt. Thus we have found and ex- 
plained the most prominent features of the stroboscopic 
many-body dynamics in our discussion of the doublet. 

Furthermore, the motion of whole clusters of fermions 
through the lattice and the exchange of fermions between 
clusters can be observed in the stroboscopic dynamics, as 
shown in Fig. |5] As expected, clusters are very stable for 
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Fig. 5. (Color online) Density plot for the time evolution of 
an initial state with clusters of different numbers of fermions 
and observation time interval J At — 2. For large interaction 
strength [panels (b) and (c)], we find clusters moving as a 
whole (indicated by triangles). Ovals indicate processes where 
single fermions are exchanged between clusters or attached to 
a new cluster. 



high interaction strengths. The hopping amplitude for 
a cluster of n fermions is of order / \ V\^~^ ^ decreas- 
ing strongly for larger clusters, as can be perceived in 
Fig-lgc). 

Expectations for other models. The previously dis- 
cussed nonmonotonic decay of a cluster is not unique 
for our model ([T]) , but can also be observed in the Fermi- 
Hubbard (FH) and Bose-Hubbard (BH) models. In the 
following we will consider clusters of doubly occupied 
sites in ID. The dissolution of a single particle from the 
edge of the cluster is well described by considering the 
survival probability of a doublon (fermionic or bosonic 
double occupancy). Using a similar analysis as for the 
doublets, one finds qualitatively the same behavior as 
shown in Fig. [4]^b). However, this description is not suf- 
ficient to understand the decay of clusters. It does not 
include the escape of paired particles, which leads to a 
complete decay of the cluster in the long time limit even 
for large onsite interaction. It also neglects tunneling 
of bosons between occupied sites. In the FH model the 
cluster is initially Pauli blocked [as in model ([T])], such 
that the dynamics is restricted to the edges. The escape 
of single fermions or doublons can be well approximated 
by the Hubbard dimer model, see [23 . The associated 
probability is ^m^{^WV{Uft/2)/[2 + 2(7/4^)^] for 
single fermions (at times Jt <C 1) and for doublons it is 
[Jtf for Jt < J/V and 4.{jH/Vf at J/V <C Jt <C 1. 
We numerically verified that this oscillating behavior of 
the single fermion escape leads to nonmonotonic decay of 
the cluster for onsite interaction strengths V/J > 4, see 
Fig. [6]^a). The probability for an initial bosonic cluster 
configuration of N sites to be destroyed is for short times 
given by 4[1 + 3(7V - l)](Jt)^. The first term in brack- 
ets stems from a single boson escaping from the edge, 
while the second and larger contribution comes from a 
boson tunneling between occupied sites. For larger times 
we studied numerically the destruction of clusters [see 
Fig. |6jb)] and find that the survival probability of the 
initial cluster configuration displays a nonmonotonic be- 
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Fig. 6. (Color online) Destruction of a cluster in the ID Fermi- 
Hubbard (a) and Bose-Hubbard model (b). We show the 
probability to detect a cluster of five doublons, see insets, 
in its initial state as a function of time for different onsite 
interaction strengths V . 



havior for V/ J >Q with a pronounced local minimum at 
t ^ 71 /V (about half the oscillation period for the dissolu- 
tion of a boson and tunneling of a boson between doubly 
occupied sites) , and the cluster is destroyed primarily via 
"internal" hopping processes. Both models show a strong 
nonmonotonic behavior of the cluster decay probability 
for large onsite interaction strengths. Thus, by tuning 
the time interval between observations in a stroboscopic 
measurement, one is able to enhance and suppress the 
decay of clusters, or even different decay channels of the 
cluster. 

Experimental realization and outlook. The Hamil- 
tonian ([T]) is related to the Heisenberg XXZ model 
by Wigner- Jordan transformation. The stroboscopic 
dynamics is identical for both models as the outcome 
of observations depends only on spatial density-density 
correlations. These Hamiltonians can be experimentally 
realized in optical lattices with fermionic polar molecules 
[26 or 2-species fermions or bosons in the insulating 
phase |27l [28] . For both realizations single-site detection 
has not yet been implemented, but experimental progress 
is being made toward this goal. Experimentally, the 
most challenging step needed to observe the interplay 
of many-body dynamics and measurements discussed 
here would be to make the observations nondestructive, 
whereas currently atoms are heated into higher site 
orbitals and atom pairs are lost due to light-induced 
collisions Beyond the scenarios discussed here, 
one may also be interested in the influence of external 
driving or measurements that either are weak or target 
only specific sites. 
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SUPPLEMENTAL MATERIAL 



Numerical implementation of the single-site resolved observation 

An ideal single-site resolved observation of particles is a projective measurement in the basis of many-particle 
configurations (occupation number states in real space). However, due to the exponentially large number of states, we 
need a numerically efficient way to sample such outcomes. This is achieved by generating the measurement outcomes 
in a stepwise fashion, building on the fact that (positive) n-particle densities pn factorize into conditional probabilities 

n 

Pn{si,...,Sn) =pl(si) • ]^Pi(Si|Si_i,...,5i). (5) 

i=2 

Here, Si G {!,..., A/'g} denotes the site of the zth particle and p^(s^|s^_i, . . . , 5i) is the conditional probability of 
finding the zth particle at site Si given that there are i — 1 particles at the sites si, . . . The procedure starts 

by randomly drawing the position of the first particle from a distribution given by the one-particle density. In the 
next step, we draw the position of the second particle, conditioned on the location of the first one, and continue 
iteratively. This way less than n • Ns values of joint probability densities have to be calculated, in comparison to the 
full number (^^) of fermionic many-body configurations. This approach relies on being able to calculate efficiently 
both the pure time evolution between observations and the z-particle densities (1 < ^ < n). In the present work, we 
use the time-dependent density-matrix renormalization group, which is an extremely powerful method for interacting 
one-dimensional systems. For the model considered in the main article, it is numerically even more efficient to draw 
the position of the first particle as before and then project the state onto those configurations where a particle is 
present at the selected site. After rescaling the resulting state, the new one-particle density is calculated. From this 
distribution we draw the position of the second fermion, excluding all sites already occupied by a fermion,and iterate 
the steps for the remaining fermions. For a discussion about local measurements on quantum many-body systems 
using matrix product states, see [Ij. 



Stroboscopic dynamics of noninteracting fermions 



In the noninteracting case, we can explicitly calculate the n-particle densities needed to simulate the stroboscopic 
dynamics. After each observation, the many-particle wave function is a Slater determinant of single-particle wave 
functions, and for noninteracting fermions this remains true even during the subsequent time evolution. Using Wick's 
theorem, the n-particle density of N fermions at time t can be written as (see also [^): 



Pn ("^1 ") ' • ' ") 1 ^) 



(6) 



Here, Ms^si{^) = XljLi '^skmj{t)Tgj^^.{t)^ with the fermion propagator Tki{t) = i^^~^^ Jk-i{2Jt)^ and the sum is taken 
over all initially occupied sites rrij. Note that the one-particle density pi{t) is just the sum of the densities of the 
individual fermions, cf. Fig. 1(b) in the main article. Motion in arbitrary potentials would be captured by different 
propagators Tki{t). 
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